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Abstract 

QCD predicts that the interactions between quarks and gluons change from a 
confining to a screened Yukawa form above a critical temperature Tc ~ 150 MeV. In 
this talk, I review some of the key observables in heavy ion reactions which are being 
used to search for this new partonic Yukawa phase at SPS and RHIC. These include 
collective observables such as dE±/dy(Pp±, meson interferometry, jet quenching, and 
J/ip suppression. 



* Talk to be published in the Proc. of 14th Nishinomiya- Yukawa Memorial Symposium Nov. 1999, 
Japan; updated with a critique of the CERN press release Feb. 8, 2000 
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§1. The Ubiquitous Yukawa 

In 1935 H. YukawaB) proposed a theory of nuclear forces based on the exchange of a 
massive boson 

VY{r,m) = aeff^— (1-1) 

He estimated that m ~ 100 MeV to account for the short range ~ 2 fm, and in 1947 Powell 
discovered the pion and confired this theory. An important theoretical precursor was the 
Klein-Gordon equation &, whose static Green's function is ( |1-1| ). Elaboration of Yukawa's 
meson theory since then (including spin, orbit, isospin vertex factors) forms the basis for the 
present effective theory of nuclear forces&'S^. 

In the unrelated field of electrolytes, Debye and Hiickle had also come across the Yukawa 
potential in another context. It emerged from solving the Poisson equation in a conducting 
medium 0*. Assuming that local charge density fluctuations occur with a Boltzmann proba- 
bility, exp{—q(f){x) /T), the polarizability of the medium in the presence of an external charge 
density, leads to a non-linear self consistent equation 

V0(x) = -47r(pe.(x) + E gn,e-^^(^)/^) . (1-2) 

In the linearized approximation, the solution for a point charge is again eq. { \L ■ 1|) , but in 
this case the effective mass is the Debye electric screening mass 

/z2 = 47r^gV,|/T . (1-3) 
q 

A dense conductive medium therefore transforms Coulomb into Yukawa. 

In nuclear theory, the Yukawa meson mass results from the flnite gap of the elementary 
excitations (pions, ...) of the physical QCD vacuum. By analogy, it should be possible to 
modify the nuclear Yukawa potential by increasing the nuclear density or temperature. In 
this talk I discuss current efforts to try to manipulate nuclear matter in the laboratory to 
force a breakdown of Yukawa's hadronic theory. As we will see, however, it seems difficult 
to escape from Yukawa. In the new deconflned, chirally symmetric phase of QCD at high 
temperatures, the Debye-Huckle mechanism takes over and the Yukawa potential between 
nucleons mutates into a color-electric screened Yukawa potential between partons of a quark- 
gluon plasma. 

In QCD, the color potential between partons is approximately Coulombic at small dis- 
tances due to the asymptotic freedom property of non-Abelian gauge theories. However, 
below a critical temperature, ~ 150 MeV, the conflnement property of the nonperturba- 
tive QCD vacuum allows only composite color singlet objects (hadrons) to "roam freely" in 
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the laboratory. The effective potential between the colored partons has a long range linear 
confining term nr to prevent them to roaming more than 1 fm away from any color neutral 
blob (hadrons). The restraining force or "string" tension, k ~ 1 GeV/fm, is huge. In this 
confining phase of QCD, the heavy qq potential is well parameterized by the Liischer formB 

VL(r,0) = -— + ft:r , (1-4) 
r 

(as long as dynamical quark pair production is ignored). The Coulombic part, with strength 
ai = 7r/12, arises from the zero point quantum fluctuations of the string (see p. 803 of 
for an intuitive derivation). This conflning potential between heavy quarks has been 
directly observed numerically using lattice QCD techniques^. As the temperature increases, 
but remains below the deconflnement transition, T < Tc, the enhanced fluctuations due to 
thermal agitation of the string modifles the effective potential into the approximate Gao 
formli* 



Vb(r,T) = -ai 



2 

-tan"^(2rT) 

TT 



TT r^o , 2 1,1, 

-T^ 1 - -tan"^ — - 
3 ^ vr ^2rT' 



r + ^ln(l + (2rT)2) 
(1-5) 



The decrease of the effective string tension, k(T), predicted above has been also observed 
in lattice QCD calculations 0* as shown in Fig.(l). However, the "measured" string tension 
is found to decreases faster than predicted in eq.(l-5) near the critical temperature. Note 
that in Fig.(l) (5 = 2Njg^ = (3, + log(T/T,), i.e., T ^ T,exp(ll/6(/3 - (3,)) with 

Pc = 4.0729 for this lattice calculation. 

For temperatures above Tc, we see from Fig.l that a new Yukawa phase of QCD is 
predicted, and that the heavy quark potential mutates into a short range generalized Yukawa 
form, which on the lattice is measured in the form 

Mr,T,d)^-^,e->'^^->"^ . (1.6) 

Note that di = 2, ^l{T) = 2mE{T) correspond to a pure Yukawa interaction in terms of the 
lattice QCD flt parameters {dL,fiL)- The perturbative thermal QCD chromo-electric Debye 
mass niE = fi{T)/2 is© 

m^(T) = g{T)T + (1-7) 

for Nc colors and Np flavors. For Nc = 2,Nf = in Fig. (1), we expect ^il{T) = 2mE = 
1.6g{T)T as shown by the solid line in Fig.(2). The flts to the lattice QCD measurements 
Fig.(l b) from0) show that for T > 2Tc jji ~ 2.5T is not far from the pQCD estimate. 
However, the exponent dL ~ 1-5 is signiflcantly below the value 2 expected from pQCD. 
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Fig. 1. (a) The heavy quark potential in the confined phase of SU(2) quenched QCD (left side) 
compared to that in the deconfined Yukawa phase above Tc (right side). Results are for a 32^ x 4 
lattice from Karsch et al0. Eq.(l-5) fits the confined lattice "data" well, but the QCD string 
tension decreases more rapidly near Tc- For T > the potential is screened by the deconfined 
gluons (in this quenched calculation) and acquires the generalized Yukawa form ( |l-6[ ). Here 
T ^ Tcexp(ll/6(/3 - (3c)) with (3c = 4.0729 and the lattice spacing is a = l/AT. 

An even more striking nonperturbative deviation is seen near Tc, at which point d < 1 and 
H ~ T/2. This suggests a rather long range interaction that may be the precursor of the 
confinement transition. This strong deviation from the perturbative QCD near Tc was also 




Fig. 2. The chromo-electric Debye screening mass, /u(T) = 2mE{T) and the effective exponent, 
d{T), of the effective Yukawa potential in the deconfined phase versus T/Tc is shown from 
Karsch et alS for SU(2) quenched QCD. The lines show expected dependence based on thermal 
pQCD. Note that near Tc, the range 2/ is much larger than predicted by pQCD and that d ~ 1 
implies an especially long range interaction there. 

reported previously in ref.0), where in addition it was found that the effective coupling is 
rather small a ~ 0.15 above Tc. It is important to keep in mind that the above numerical 
experiments do not include dynamical quarks and are limited to SU(2). Nevertheless, they 
provide strong evidence that QCD predicts a qualitatively new (nonperturbative) partonic 
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Yukawa phase of matter that should exist at an energy density only an order of magnitude 
above that in ground state nuclei (e > 2 GeV/fm^). 

The thermodynamic properties of the deconfined QCD phase are shown in Fig. 3 from 
ref.lli' for 2 flavor 12^x6 IQCD. A present limitation of all lattice results so far is that the pion 
is still too massive to make make contact with the "known" thermodynamic properties of 
ordinary hadronic/nuclear matter below T^. Recent advances in implementing Domain Wall 
Fermions on the lattice and the availability of new TeraFlop scale computers at Columbia and 
the Riken Brookhaven Research CenteM and the CP-PACS project in Tsukubalii', should 
enable much more precise calculations of the quark-gluon plasma equation of state in the near 
future. Two striking features of Fig. (3) suggest two key observable signatures of this phase 



Fit using 0(4) scaling 




T(GeV) T(GeV) 

Fig. 3. Thermodynamic energy density (e/T^ top curves left), pressure {Sp/T"^ lower curves left) 
and speed of sound squared (right) from lattice QCD (2 flavor 12^ x 6) from the MILC collablli'. 
The curves are zero quark mass extrapolations. Note the rapid reduction of the pressure and 
speed of sound as Tc is approached from above. 



transition in nuclear collisions. First, the entropy density o"(T) = {e + p) /T increase very 
rapidly with T in a narrow interval AT/T^ < 0.1. Second, the plasma becomes extremely 
soft p/e <^ 1 and ^ 1 near Tc. As we review below, the first feature can lead to time 
delay (the QGP stall) measurable via hadronic interferometry. The second feature can lead 
to interesting non-linear collective flow observables in nuclear collisions. The experimental 
verification of these fundamental predictions of QCD is the primary motivatiion for the heavy 
ion experimental program at Brookhaven and CERN. In the following sections, I review first 
how deep into the QGP phase RHIC may be able to reach, and then discuss several signatures 
used to test the QCD predictions in such experiments. 
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§2. Initial conditions in A+A 

In order to see the partonic Yukawa phase, we must first create an extended blob of matter 
at 100 times the density of nuclear matter. Head on collisions of heavy ions are used for that. 
There are many ways to describe how this dense matter is formed. One intuitive picture 
is given in terms of the McLerran-Venugopalan models'*. For highly boosted nuclei with 
Ecm ~ 100m AT, time dilation effectively freezes out the quantum chromo fiuctuations inside 
the nuclei while the two pass through each other. Heavy Au beams can then be regarded 
as well collimated, ultra dense beams of partons. This (chromo Weizsacker- Williams) gluon 
cloud contains a very large number, Ga{x,Pq) ~ A/x^^^, of almost on-shell coUinear gluons 
with longitudinal momentum fraction x = po/Ecm <^ 1- As the clouds pass through each 
other, partons scatter via chromo Rutherford and decohere into a mostly gluon plasma on a 
fast time scale l/p± ^ 1 fm/c. The number of gluons pairs (mini-jets) extracted from the 
nuclei by this mechanism at rapidities yi and transverse momentum ±k± can be calculated 
in pQCD as follows 

^^^/Ig = Kx.GAxi, kDx.Gsix,, kl)^^§p^TMb) , (2-1) 
dyidy2dkj_ dkj_ 

where xi = x±{exp{yi) + exp(?/2)) and X2 = x±{exp{—yi) + exp(— ?/2)), with x± = k_}^/ ^s, 
and where the pQCD gg — > gg cross section for scattering with t = —k'^{l + exp(y2 — yi)) 
and y2 — yi = y is given by 

da^s ^ 9 47ra2 {1 + 6^ + e"^)^ 

~dr ~ 8 ki (e2^/2 + e-2^/2)6 • ^^'^^ 

The nuclear baryon number. A, only plays here the role of increasing the density of partons 
and providing the geometrical amplification, T^Bib) ~ 30/mb, for the number of binary 
nucleon-nucleon collisions that occur per unit area. A factor K ^ 2 simulates next to 
leading order corrections. 

For symmetric systems, A + A, with Ga ~ AG. the inclusive gluon jet production 
cross section is obtained by integrating over y2 with yi = y and k^ fixed. To about 50% 
accuracy, the single inclusive gluon rapidity density in central collisions can be estimated by 
the following simple pocket formulai^-* 

dN 42 /J el\ 

;2N,ix^,t)x^Gix^, t) -f ^ (2-3) 



dydt t:B? \ dt 



R 



where Ng(x±, t) = /J^ dxG{x, t) is the total number of "hard" gluons coming down the beam 
pipe in an nuclear area nR^ that can knock out unsuspecting gluons from the other nucleus. 
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This copious mini-jet mechanism is beheved to the dominant source of the gluon plasma 
that will be created when RHIC (finally) begins operation. 

Recent upper bound estimates of the total gluon rapidity density as a function of the 
CM energy from^ are shown in Fig. (4). The differential yields are integrated down to a 
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Fig. 4. Expected copious mini-jet production in A-l-A collisions from ref. 0). The dashed curves 
in both figures correspond to an estimated maximum number (saturation) at which the nuclear 
area if filled with closed packed mini jet gluons {N{po) = PqR^). {Ni is dimensionless, [mb] 
is a typo in fig.). The magnitude of possible hydrodynamic transverse energy loss due to 
longitudinal work in an ideal p = e/3 quark gluon plasma is shown by Ej^j- 

transverse momentum scale po ~ 1 ~ 2 GeV. This scale separates the "soft" nonperturbative 
beam jet fragmentation domain from the calculable perturbative one above. The curve 
marked saturationlll* is an upper bound marking the point where the transverse gluon density 
of mini-jets becomes so high that the newly liberated gluons completely fill the nuclear area, 
i.e., dN/dy ~ PqR^- At that point, higher order gluon absorption may limit the further 
increase of the gluon number. At RHIC energies these estimates yield up 1500 gluons per 
unit rapidity. My more conservative estimates together with X.N. Wangilt'S* gives 500 
gluons per unit rapidity when initial and final state radiation is also taken into account. 
This is obtained with a fixed po = 2 GeV, that was found to be consistent with all available 
pp and low energy AA data using the HIJING mini-jet event. 

The energy dependence of the final charged particle radidity density from HIJING, in- 
cluding soft beam jet fragmentation, is shown in Fig. (5). Approximately one half of the 
height "Mt. RHIC" comes from soft beam jet fragmentation processes modeled in HIJING 
using Lund/Fritiof strings. The initial energy density reached in such collisions can be 
estimated using the Bjorken formula 

1 dEr 



iro) 



(2-4) 



ttR^Tq dy 

For po ~ 1 - 2 GeV, dEr/dy ~ 400 - 2000 GeV, and so e(ro ~ 0.5 fm/c) > 10 GeV/fm^ 
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Fig. 5. Prediction of the height of "Mt. RHIC" at E^m = 100 AGeV based on HIJINgEiI' event 
generator from Quark Matter 95©. Comparison to measured charged particle rapidity densi- 
ties in Au + Au reactions from Eiab = 1, 10, 200 AGeV is also shown. 

should be easily reached, well inside the the deconfinement phase of QCD. At SPS energies, 
on the other hand, my estimates indicate that nuclear collisions may just reach the transition 
region and depart from that region in a very short time. 

Since the the time of this talk, provocative and somewhat overstated press releases^ have 
been issued from CERN stating that "the experiments on CERN's Heavy Ion programme 
presented compelling evidence for the existence of a new state of matter in which quarks, 
instead of being bound up into more complex particles such as protons and neutrons, are 
liberated to roam freely. ... We now have evidence of a new state of matter where quarks 
and gluons are not confined." 

As discussed below, I disagree with the above interpretation of that truly impressive body 
of data. What is compelling is that some form of matter, much denser than ever studied 
before, was created. Inferences about quark and gluon degrees of freedom are based on 
qualitative scenarios and schematic models. At the relatively low momentum scales accessible 
at SPS energies, the quarks and gluon degrees of freedom in the dense matter are mostly not 
resolvable. Even at the highest transverse momentum, the pion spectra were shown to be 
very sensitive to nonperturbative model assumptions such as to the magnitude of intrinsic 
momenta and soft initial state interactions The dynamics of the non-perturbative beam 
jet fragmentation and hadronic final state interactions cannot be disentangled at SFsS-*. It is 
useful to recall that in experiments on pp scattering, it was only possibly to see unambiguous 
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evidence for the tell-tale Rutherford scattering of point-like partons when collider energies 
^/s > 200 - 2000 GeV became available. At RHIC high p_L > 10 GeV probes become 
kinematically available and hence very small wavelength resolution of partonic degrees of 
freedom finally becomes kinematically possible. While there is an abundance of interesting 
signatures showing that dense matter was formed at the SPS (through the non-linear in 
dependence of several observables on multiplicity or A), the bottom line is that those data 
have said nothing about whether the QCD predictions in Figs 1-3 are correct or not. We 
simply need higher resolution. RHIC, with its factor of ten increase in CM. energy reaches 
a factor of ten deeper into the new phase reaches. Coupled with the availability of ten times 
shorter wavelength probes, it should finally become possible to actually see direct evidence 
of "freely roaming quarks and gluons" . 

2.1. Global Signatures of Collective Dynamics 

The simplest but often ignored global barometer of collectivity in nuclear reactions is 
the A and energy dependence of the transverse energy and charged particle rapidity density. 
At RHIC energies, HIJING predicts that initial transverse energy density in central A -|- A 
collisions scales nonlinearly with A 

i^«lGeV^"(l + log^) (2.5) 

This leads to about 1 TeV per unit rapidity, on the other hand, it was found that 

gluon saturation could limit the A dependence to approximately linear, E_i_ ~ y4^ °^, but with 
a value several times that of HIJING. In Fig. (4), the initial gluon density actually grows less 
than linear A"'^^ in that model. The A^'^ scaling of HIJING with its conservative po = 2 
GeV scale fixed by pp data is simply due to the number of binary interactions via ( p-3|) . The 
initial energy density in HIJING varies approximately as 

eo ^ 0.6 GeV/fm' A^'^^ ^ ^) _ ^2-6) 

In0^ its magnitude and scaling are predicted to go as e(r = 1/psat) ~ O.lA^-^s^'^^. 

If local equilibrium is achieved and maintained, then a very basic prediction of hydro- 
dynamics is that longitudinal boost invariant expansion together with pdV work done by 
pushing matter down the beam pipe will cool the plasma and convert some its random trans- 
verse energy into collective longitudinal kinetic energy. For an equation of state, p = c^e, 
this cooling and expansion causes the energy density to decrease with proper time as 

<r) = e(ro)(^)'^'' (2-7) 



9 



Entropy conservation leads to a conservation of rn^r), where n is the proper parton density. 
At least if the matter is initially deep in the plasma phase, then (as seen in Fig. 3) longitudinal 
will be done with ^ 1/3. Consequently, the transverse energy per particle should decrease 
by a factor 2-3 before freeze-out as0'''0) 

^.(r) = '-^ = e,in)r^Y- (2-8) 



dN ' \ t 

However, dissipative effects due to finite mean free paths reduce the effective pressure in any 
system. For the Bjorken expansion, the relaxation time, = l/{aTn) oc r, then increases 
with time as n decreases. Numerical solution of 3+1 D transport equations with pQCD cross 
sections, ctt ~ 2 mb, indicate that dissipation reduces the transverse energy loss for HIJING 
initial conditions rather significantly (see detailed comparisons inll^). Ideal hydrodynamics 
predicts that about one half of the initial produced transverse energy goes into longitudinal 
work 

It is important to keep in mind that the commonly assumed freeze-out prescription with 
Tf fixed on a fixed freeze-out energy density hypersurface, eijf) = e/ ~ 0.15 GeV/fm^ is 
only a rough prescription that is never accurate^. As the interaction size decreases {A 
or multiplicity decreases), the effective speed of sound due to dissipation decreases and the 
system freezes out earlier. A detailed study of the A or multiplicity dependence of dE±/dy 
and ct is needed to calibrate the interplay of pdV work and dissipation on these barometers. 

One of the major experimental discoveries of WA98, NA49 and the other SPS experiments 
is that at those energies dE±/dy as well as dN/dy scale with A or wounded nucleon number 
nearly linearly il* as shown in Fig. (6). Both the E± and the charge multiplicity increase as 
~ ^i 07_ 'jj-^ggg findings differ from the VENUS model, which has considerable nonlinearly 
due to the assumed sea string contributions. In fact, simple Glauber wounded nucleon 
models reproduce very well the nearly linear correlation between E± and the veto calorimeter 
(spectator) energy observed in all experiments at spsii). The implications of these data 
depend on the A dependence of the initial conditions of course. One view0 is that the 
initial can scale arbitrarily with A, but because perfect local equilibrium is maintained 
up to a critical sharp freeze-out hypersurface, the final E^ and always scale linearly with 
A. My view is that at the SPS, pQCD is mostly inapplicable to bulk phenomena and the 
linear dependence arises from additive nature of soft beam fragmentation together with the 
absence of pdV work at early times. If the QGP transition region is just barely reached, as 
I believe without experimental proof, then the softness of the QCD equation of state with 
cj. <^ 1/3 seen in Fig. (3) and dissipation can conspire to prevent the dense matter from 
performing longitudinal work. However, it is impossible to tell from the data whether the 
observed null effect is then due to a low pressure Hagedorn resonance gas of hadrons or to 
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Fig. 6. Scaling of transverse energy per charged particle and and the charged particle radidity 
density with number of participants in Ph + Ph from H. Schlagheck et al WA98@*. The nearly 
linear scaling in the number of participants can be interpreted as no longitudinal work was 
done by the lazy dense matter at the SPS. 



a low pressure lazy "plasma" with ^ 1. One has to go to higher energies to the plasma a 
chance to work. 

In spite of the above remarks, ideal Euler hydrodynamic calculations have in fact been 
able to fit the null effect in the data quite well. This is possible by assuming arbitrary initial 
conditions chosen such that that after hydro does its work, the final results just happens to 
reproduce the data. An illustration of this degree of freedom is shown in Fig. (7) fromS). By 
adjusting the initial four velocity field appropriately, the calculated final pion distributions 
can be made to reproduce the data starting from two completely different initial energy 
density profiles. The same good fit to data can be obtained also for 74A/^F assumed equation 
of state of dense matter simply by readjusting the initial conditions appropriately. Therefore, 
it is not surprising that even hand calculator fireball models are able to fit much of the data 
without invoking any dynamical or thermodynamical assumptions about the properties of 
dense matter prior to freeze-out 11* '0''. It is clear that inferring the existence of freely roaming 
quarks and gluons based on such fits is impossible 

In order to use the global barometers to search for evidence of collective phenomena in 
dense matter at RHIC, it is first necessary to eliminate the freedom to choose arbitrary initial 
conditions H-* . At collider energies, and non-Abelian field techniques 0^'© 

provide the needed theoretical calibration tools to fix initial conditions At RHIC and LHC, 
the initial plasma is expected in any case to be so deep in the deconfined phase that the 
soft p <^ e/3 transition region and dissipation cannot spoil the longitudinal work as is the 
case at SPS. From detailed covariant transport calculations i^-^'S* the barometric evidence 
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for longitudinal work should finally be observable in spite of finite mean free path effects. 
Deconvolution of the equation of state and dissipative corrections requires however a detailed 
study of the A and multiplicity dependence of the global Et barometer. 




Fig. 7. Some different initial conditions that lead with suitable initial velocity fields via ideal 
hydrodynamics to same final observable pion rapidity density at SPS from ref.iil). 

This is the first important signatures to look forward to at RHIC. While it is not possible 
to converge on why no longitudinal work was observed at the SPS, a positive signature at 
RHIC would render such discussion mute. We can only be sure that a plasma was created 
if it does something collective! The global Et and multiplicity systematics as a function of 
centrality as well as A therefore provide key handles in this search. 



§3. Transverse Flow 

In contrast to global transverse energy barometer a completely different measure of baro- 
metric collectivity is afforded by the study of the triple differential distributions, dN/dy(Ppj_. 
Already at sub- luminal Bevalac energies (< 1 AGeV), azimuthally asymmetric collective di- 
rected and elliptic flow were discovered long ago. For non central collisions, 6 7^ 0, the 
asymmetric transverse coordinate profile of the reaction region leads to different gradients 
of the pressure as a function of the azimuthal angle relative to beam axis. This leads to 
a "bounce" off of projectile and target fragments in the reaction plane and to azimuthally 
asymmetric transverse momentum dependence of particles with short mean free at mid ra- 
pidities. This phenomenon has now been observed at both AGS and SPS energies as well. It 
will certainly be there also at RHIC and LHC. In Fig. (8) the first two Fourier components 
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of the azimuthal flow patterns are shown: 
dN 



vo{l + 2vi cos(0 - (j)R) + 2v2 cos(2(0 - 0^)) + ■ ■ ■) 



(3-1) 



dyd'^p± 

Azimuthal asymmetric collectivity is clearly observed at the SPS. The important question 
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Fig. 8. Qeft) Azimuthal flow analysis in Pb + Pb from H. Schlagheck et al WAOsi!''!!) and NA49 
datacj). (Right) Predicted dependence of V2 at mid rapidity Au + Au at RHIC from parton 
cascade calculations ZPC 

i) with HIJING initial conditions. 



is how this type of barometer could serve to help search for evidence of the QCD transition. 
Unlike the global E± barometer discussed in the previous section, transverse flow can develop 
at later times because the gradients are controlled by the transverse size of the nucleus, not 
the proper time interval relative to formation. InS-' the idea was proposed that one could 
use V2 for example to study the predicted softening of the QCD equation of state. Typically, 
hydrodynamics calculations lead to a factor of two smaller V2 for an equation of state with 
a soft critical point as in Fig 3 versus one in which the speed of sound remains l/Vs. 
Searches for anomalous V2 dependence as well as vi are underwayElt. As with the global 
barometer, dissipation can of course also simulate a soft equation of state. In ref@* we 
studied the dependence of V2 on the transport parton cross section. The results shown in 
Fig(8) for HIJING initial conditions indicate that there is indeed a significant reduction of 
V2 relative to hydrodynamics, but that the impact parameter (multiplicity) dependence of 
that observable can again be used to disentangle the equation of state versus viscosity effects 
as in the case of the global Et barometer. 
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4. The QGP Stall and Time delay 



Hadron interferometry has been developed into a fine art in heavy ion colhsions to image 
the space-time region of the decoupling 4- volume. In@* it was proposed that a possible 
signature of the QGP phase transition would be a time delay associated with very slow 
hadronization. The plasma "burns" into hadronic ashes along deflagration front that moves 
very slowly if the entropy drop across the transition is large©. Fig.9 shows the evolution of a 
Bjorken cylinder with time and transverse coordinate from a hydrodynamic simulation with 
different equations of state from©. The main point to note is that time delay is a robust 
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Fig. 9. Evolution of isothermal contours for a boost invariant QGP cylinder with initial energy 

density ~ 20 GeV/fm'^ for different equations of state from ref.@*. For a strong first order 
transition the QGP "stall" is very strong, but it also remains significant even if the transition 
is only a rapid cross over consistent with Fig. 3. For an ideal gas with p = e/3 on the other 
hand, the freeze-out time scale and radius are approximately equal. An experimental probe of 
this decoupling geometry is provided by pion or kaon interferometry. 

generic signature of a rapid cross over transition of the entropy density. In particular the 
"stall" is expected even for a smooth cross-over transition as long as the width AT /T^ < 0.1. 
However, its magnitude also depends on the entropy drop across that region. The figures are 
for an entropy drop by a factor 10 consistent with the lattice QCD results. Unfortunately, as 
noted before lattice QCD has not yet resolved the hadronic world below due to numerical 
problems. If the entropy jump is much smaller, then this signature would also disappear. 

High statistics measurements of pion and kaon interferometry searches for time delay at 
AGS and SPS have come up empty handed thus far. No time delay has ever been observed in 
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any nuclear reactions thus far. This could be due to (a) the absence of a large rapid entropy 
drop in real QCD, or (b) to unfavorable kinematic conditions at AGS and SPS energies. From 
the hydrodynamic calculations it was found that a large time delay signal requires that 
the initial energy density be deep within the plasma phase. The slow deflagration waves seen 
in Fig. 9 do not arise if the initial energy density is close to the transition region because the 
initial longitudinal expansion cools the plasma too rapidly. The optimal conditions to see 
this effect was predicted© to occur at RHIC energies. For much higher energies (LHC), on 
the other hand, the transverse expansion of the plasma has too much time to develop and 
that spoils the possibility of a slowly burning plasma stall. I note that further work© has 
shown that high kaon interferometry adds an especially sensitivity handle in the search 
for this time delay signature. If ever observed, the time delay signature would be smoking 
gun that a new form of matter with bulk collective properties was created. 



§5. The J/Psi puzzle 

In 1986, Matsui and Satz proposed an intriguing direct measure of the transmutation of 
the qq forces in Fig. 1. The idea was that J/tp can form in the vacuum because the confining 
Luscher potential can bind a cc pair into that vector meson. If that pair were placed in 
a hot medium in which the chromo-Debye screening potential is short ranged, then above 
the temperature where the screening length is smaller than the J/tp radius, the cc would 
become unbound and and the charm quarks would emerge from the reaction region as an 
open charm DD pair. They thus predicted that J ftp suppression would be a smoking gun 
for the deconfinement transition. 

J/ip suppression was first seen in 1987 in O + U reactions by NA38. Since then this 
smoking gun has (unfortunately) never stopped smoking! J/ip suppression seems to be as 
ubiquitous as the Yukawa potential. It is now clearly observed in p + A reactions as seen in 
Fig. (10). High mass Drell-Yan pairs, on the other hand, formed via qq ii was observed 
to scale perfectly linearly with the number of binary collisions. This is because lepton pairs 
suffer no final state interactions and the quark initial state (Cronin) interactions are invisible 
in p±_ integrated DY yields. 

J/ip are suppressed on the other hand as (AB)^'^ due to some nuclear effect that is 
independent of QGP production. In the recent Pb+Pb analysis an excess 25% suppression 
of J/ip was observed. NA50 claims that this enhanced suppression relative to the (AB)^'^ 
trend from lighter projectile p,0, S + U data is finally the real smoking gun0\ From the 
suppression pattern in Pb as observed as a function of centrality (see Fig. (11)), NA50 has 
now claimed that in fact "Together with the results previously established by the NA38 and 
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Fig. 10. Drell-Yan scales as {ABy very accurately from NA38, NA50. However, J/^p scale as 
(AB)^'^ showing that their production is suppressed with increasing nuclear volume even in 
p + A reactions Ref. The "anomalous" suppression in Pb + Pb is the 25% deviation from 
the empirical (AB)^'^ expectation. 
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Fig. 11. "Evidence for deconfinement of quarks and gluons from the J/psi sup pr ession pattern 
measured in Pb Pb collisions at the CERN-SPS" claimed by NA50 in Ref.cS^. The curves 
on left show transport theory estimates for hadronic final state dissociation. The enhanced 
suppression relative to exp{—aNPoL) nuclear suppression observed inp + A and S + U is shown 
versus a rough estimate of the initial energy density (this scale is uncertain to at least a factor 
of two for each point!). 



NA50 collaborations, a rather clear picture emerges, indicating a step-wise pattern, with no 
visible saturation in the collisions generating the highest energy densities and temperatures. 
Our observations exclude the presently available models of J ftp suppression based on the 
absorption of the J/ip mesons by interactions with the surrounding hadronic (confined) 
matter. The first anomalous step can be understood as due to the disappearance of the 
X mesons, responsible for a fraction of the observed J/ ip yield through its (experimentally 
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unidentified) radiative decay. In proton induced collisions this fraction is around 30-40%. 
The second drop signals the presence of energy densities high enough to also dissolve the 
more tightly bound J /if) charmonium state." 

While the deviation from the empirical (y4i?)°'^ scaling is very clear, the dynamic origin of 
the effect is far from clear in my opinion. As shown by the several curves in Fig. (11), a ~50% 
drop in the J/ yield as a function of Et is consistent with final state co-moving hadronic 
absorption. While the details are not reproduced accurately, large theoretical uncertainties 
about several key dynamical ingredients preclude precise comparisons at this time. It is 
important to emphasize that the so called "conventional hadronic" models suffer from just 
as large theoretical uncertainties as the plasma scenario models. Key uncertain elements 
include (1) the dynamical treatment of cold nuclear absorption responsible for the AB^'^ 
suppression even in p + A, (2) the unknown hadronic M + DDX reaction rates, and 
(3) the actual density evolution, e(x, t) needed in order to make more precise calculations. 
Only a schematic "octet model" of pre-hadronization cc interactions in cold nuclei has been 
used in the above transport models to address the first effect. That this is highly uncertain 
is shown in the work of Ref.H). The observed J/ip production cross section even without 
final state interactions is suppressed as follows: 



where F is the formation probability of the ip from a cc that emerges from the cold nuclear 
target with an invariant mass < 4M|,. If the pre-resonance cc pair multiply scatters in the 
nucleus random walk would increase + 5q^{apL) linearly with nuclear thickness as 

shown by the curve G in Fig. (12). This leads to an approximate exponential suppression that 
can be fit well by an approximate Glauber nuclear absorption factor ansatz, exp{—aeffpL) 
if a Gaussian assumption is made. This Gaussian model can thus account for the observed 
{AB)^-^ scaling light projectiles. However, it was shown inS that if the power law tails 
due to induced radiation in the medium is included (resulting from the multiple Rutherford 
rescattering of the color octet cc) , then an additional nonlinear suppression in the nuclear 
thickness L could result (curve P). This is because radiation provides another way to increase 
the invariant mass of the pre-resonance cc that can further reduce the probability for the pair 
to fit inside the ip wavefunction. While this model is also schematic, I feel that it captures 
an important element of pre-formation transient dynamical effects in nuclei that can serve 
as additional sources of nonlinear suppression, without even considering the sought after 
suppression in the comoving dense matter. The transport model curves in Fig. 11 co-moving 
dissociation possesses and pre-resonance nuclear dynamics as in Fig. 12 should be combined 
in future analysis searching for the dynamical origin of beyond nuclear Glauber absorption 
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Fig. 12. Sensitivity of "pre-resonance" J/V' nuclear absorption to details of color octet models 
of cc interactions in a nuclear medium from Qiu et altlP. The curve marked P for power law 
accounts for the anomalous absorption seen in Pb+Ph and deviated from Gaussian Glauber-like 
G expectations. 

in Ph + Ph. 

Therefore, the claim of anomalous suppression cannot rest therefore merely on generic 
enhanced suppression in Ph. It must and has been made on the possible existence of singular 
"step-like" structure of the suppression pattern. The evidence for "steps" is however the 
weakest link experimentally because a rigorous test including the substantial systematic 
errors in the Et scale has yet to be performed. This very difficult experiment has fought 
valiantly for years to reduce the systematic errors associated with calibration of Et, thick 
target multiple interactions, and calibrating the Drell-Yan and other backgrounds. This is 
illustrated in Fig. (13) where plasma scenario to 1997 data are also shown. In Fig. (14) 

recent plasma scenario fitsii to the latest data show how much the data has evolved (see© 
for details). It is also clear from Figs. (13, 14) that the plasma scenario step-function models 
while capable of fitting the data cannot be be considered as prediction of QCD. Percolation 
models serve to motivate the steps. 

Nevertheless, if the "step-like" suppression pattern survives further experimental scrutiny, 
it would certainly be the most dramatic nonlinearity observed at SPS. In this connection, 
it is also important to remark that the energy density scale in Fig. (11 b) is uncertain to 
at least a factor of two. There is no direct experimental measurement of the initial energy 
density. The energy scale there is infered from a RQMD model calculation. However, as 
shown in Fig. (7) a nice hydro fit to the data could be obtained with an initial energy density 
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Fig. 14. Percolation model step-function fits to the survival propability normalized to the em- 
pirical nuclear suppression effect from H. Satz, Quark matter 99©. Open square NA50 data 
are from thinner target runs with Pb + Pb. 



ten times higher than assumed in Fig. 11 b. 

Given the boldness of the NA50 claims, it is important to scrutinize both the experimental 
and theoretical foundations on which those claims are made. In science guilt is to assumed 
until innocence is proven. Based on the previous discussion, I remain skeptical. 

An additional problem with the plasma scenario interpretations can be seen form the 
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observed Et dependence of the J/ip Pl spectra in Fig.(??). Standard Glauber multiple 
collisions lead to a random walk in transverse momentum that are expected to lead to0^ 

{pI)ab = {pI)pp + j5pI (5-2) 

his is found to hold^-* in all reaction including Pb + Pb. In contrast, in the plasma sce- 
nario 

i) 

, only those are expected to survive that are near the surface where the nuclear 
depth L is small. Thus the prediction as shown in Fig(??) was that the should begin 
to DECREASE with increasing Ex. This was not observed^*'©. Clearly, much more work 




Fig. 15. The mean transverse momentum of the J/'0e3) shows clear evidence of multiple scattering 
in the nuclear target as expected from Glauber theory. Even in Pb + Pb the increase of the 
Pt is understood from the same mechanism in contradiction^ to predictions based on the 
plasma scenario, where at high Et, the surviving -0 are expected to have been produced only 
near the nuclear surface regions i). 



is required to sort out the very interesting suppression pattern observed by NA38,NA50. 
Experimentally, the claims would carry considerable more punch if similar " step- wise" pat- 
ternswere observed in other systems ,e.g. Xe + Xe suitably shifted in Et due to the expected 
smaller energy densities achieved. While it is not likely that such further measurements will 
get done, they will be at RHIC. A clear prediction by H. Satz at QMggii, was that under 
RHIC conditions of higher energy density, the same step wise pattern as in Fig. (14) should be 
observable in Cu + Cu interactions. PHENIX will provide a definitive test of this prediction 
soon. 
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6. The High p±^ Frontier 



One of the new areas that RHIC will open in the experimental search for the next Yukawa 
phase of QCD is high transverse momentum (short wavelength) jet probes. The rates of jet 
production and its fragmentation in the vacuum are well understood. The new physics here 
is the study of partonic interactions at extreme densities through the phenomenon of jet 
quenching S-* . Final state interactions of a jet in a dense QGP are expected to induce a large 
radiative energy lossill'. In fact, it was discovered in by BDMPsil'' that non-Abelian energy 
loss is in fact non-linear as a function of the thickness of the medium. Tests of this and other 
aspects of non-Abelian multiple collision dynamical will be possible at RHIcIl^. 

At SPS energies, this physics is out of reach given the dominance of nonperturbative 
effects as shown in Fig. (16) from ref.ii^. HIJING accidently fits the WA98 data with or 
without jet quenching. At the SPS no clean separation of soft and hard dynamics is kine- 
matically possible. However, at RHIC energies, the power law tails of the single inclusive 
distributions stick far enough above the confusing soft "noise" to gain sensitivity to the 
form of the non-Abelian dE/dx as seen in Fig. (16). This problem is closely related also 
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Fig. 16. Jet quenching at SPS vs RHIC from ref.ill compared to WA98 data. At RHIC the power 
law tail extends finally far enough above the nonperturbative "noise" to make jet quenching 
observable. 




to the problem of pre-resonance J/ip absorption discussed previously Understanding 
jet quenching is also prerequisite in testing the dynamical assumptions of recent covariant 
parton transport theories il-^'il*. 
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7. Summary 



The next Yukawa phase is awaiting discovery. The SPS data have provided many in- 
triguing indirect hints that new physics is operating in dense matter. Many puzzles, claims 
and counter claims remain because at SPS energies, both hadronic and partonic models 
have partial overlapping domains of vahdity. This "duality" is analogous to the problem 
of interpreting the R factor in e^e" coUisions below -y/s < 10 GeV. The ratio of hadronic 
to leptonic production cross sections only reaches the the magic 11/3 of pQCD above that 
threshold region where the ss, cc, bbT vector mesons dominate the nonperturbative hadronic 
physics. Similarly SPS is at the door-step where hardonic resonances begin to melt away 
and a pQCD continuum description start to become more relevant. 

With RHIC the factor of ten increase in the initial energy density will be unambigously 
in the QGP continuum. The matter so formed will also have much more time to develop 
collective signatures. The factor of ten smaller wavelength probes will finally allow experi- 
mentalists to resolve (i.e. see) the quark and gluon degrees of freedom of that plasma. Direct 
observation of longitudinal work, transverse azimuthal collectivity, time delay, step- wise J/ 
suppression in Cu -\- Cu, and jet quenching will offer direct signatures of the sought after 
new phase of QCD matter. 
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